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We study inflation and preheating in F(7Z) supergravity characterized by two mass scales of a 
scalar degree of freedom (scalaron): M (associated with the inflationary era) and m (associated with 
the preheating era). The allowed values of the masses M and m are derived from the amplitude of the 
CMB temperature anisotropics. We show that our model is consistent with the joint observational 
constraints of WMAP and other measurements in the regime where a sufficient amount of inflation 
(with the number of e- foldings larger than 50) is realized. In the low-energy regime relevant to 
preheating, we derive the effective scalar potential in the presence of a pseudo-scalar field \ coupled 
to the inflaton (scalaron) field (f>. If m is much larger than M, we find that there exists the preheating 
stage in which the field perturbations 8\ and 8(f) rapidly grow by a broad parametric resonance. 
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I. INTRODUCTION 

The gravitational theory dubbed f{R) gravity, where 
the Lagrangian / is a function of the Ricci scalar R, is the 
viable phenomenological framework to describe inflation 
and reheating in the early Universe [l[ (see also Refs. 0, 
jH for early works). For the simple model described by the 
Lagrangian f(R) = R—R 2 /(6M 2 ), where M corresponds 
to the mass of a scalar degree of freedom ("scalaron"), 
the presence of the quadratic term R 2 /(6M 2 ) leads to 
cosmic acceleration followed by a gravitational reheating 
0] • Moreover the predicted power spectra of scalar and 
tensor perturbations in this model are compatible with 
the observations of the CMB temperature anisotropics 
(see Refs. [(HI] for reviews). 

An embedding (or derivation) of a viable f(R) grav- 
ity model from a more fundamental physical theory nat- 
urally leads to supergravity (the theory of local super- 
symmetry) as the first step, because supersymmetry is 
the leading and well motivated proposal for new physics 
beyond the Standard Model of elementary particles. Su- 
pergravity is also the low-energy effective field theory of 
superstrings. 

A manifestly Af = 1 locally supersymmetric (minimal) 
extension of f(R) gravity was recently constructed in a 
curved superspace @. This was dubbed F(7t) super- 
gravity, where F(1Z) is a holomorphic function of the 
covariantly-chiral scalar curvature supcrficld 1Z. The ba- 
sic features of F(1Z) supergravity and some of its non- 
trivial models were systematically studied in Refs. (Tol — 

The first phcnomenologically viable inflationary model 
based on F(1Z) supergravity was proposed in Ref. [ijj. In 
the high-energy regime relevant to inflation this is similar 
to the Starobinsky's f(R) model with the quadratic term 



-i? 2 /(6M 2 ) Q, but the correction of the form (-R) 3 / 2 
is also present. The viability of this model was proven in 
certain limit of its parameter space, but phenomenolog- 
ical and observational bounds on the parameter values 
were not found yet. 

In the low-energy regime relevant to reheating the pro- 
posed model of F(1Z) supergravity [l4| is approximately 
described by the Lagrangian f(R) = R — i? 2 /(6m 2 ), 
where the mass scale m is not identical to the mass M in 
the high-energy regime. In this case the dynamics of re- 
heating should be different from that studied in Refs. [3| 
for the original Starobinsky's f(R) model. It was conjec- 
tured in Ref. [l4[ that the model [3] gives rise to efficient 
preheating after inflation, but it was not proven or veri- 
fied by concrete calculations. 

In this paper we extend the analysis of Ref. [3] by 
providing phy sical interpretation to the parameters of 



the model LJJ , confirm its consistency and viability, and 
find observational bounds on its parameters. In the low- 
curvature regime we also derive the effective scalar po- 
tential of a scalaron field <j> coupled to a pseudo-scalar 
field x- This potential is employed for numerical analysis 
of the preheating stage after inflation. We show the ex- 
istence of a broad parametric resonance, by which both 
the field perturbations 6x and Sip are amplified in the 
parameter region where m is much larger than M. 

Our paper is organized as follows. In Sec. |TT] we de- 
scribe how f(R) gravity arises from F(1Z) supergrav- 
ity and review the model fl4j by adding more details. 
In Sec. IIIII we study the inflationary dynamics in the 
model [14j under the slow-roll approximation and pro- 
vide phenomenological bounds on the model parameters. 
In Sec. II VI our model [l4[ is confronted with the observa- 
tional tests of CMB and other measurements. In Sec. |V] 
we compute the effective scalar potential of two physical 



2 



scalars in the low-energy regime. In Sec. IVII wc numeri- 
cally study the dynamics of preheating of the multi-field 
system. Our conclusion is given in Sec. IVIII 



II. MODEL 

The action of F(1Z) supergravity in the chiral (curved) 
J\f = 1 superspace of (1 + 3)-dhncnsional spacetime is 
given by 



S 



d 4 xd 2 6£ F(TL) +H.c. 



(1) 



where F(1Z) is a holomorphic function of the covariantly- 
chiral scalar curvature superfield 1Z, and £ is the chiral 
superspace density [l5[. The scalar curvature R appears 
as the field coefficient of the 9 2 term in the superfield 
1Z. We use the metric signature (+,—,—, — ), so that the 
sign of R is opposite to that of Ref. Q . See Ref. [IH for 
more details about our notation and F(1Z) supergravity. 

The chiral superspace density (in a Wess-Zumino type 
gauge) reads 



z{x) [1 - 2ia a ^ a (x) +9 2 B(x)] 



(2) 



where e = y/— det g^, g^ w is a spacetime metric, ip% = 
e^ipQ is a chiral gravitino, and B = S — iP is the complex 
scalar auxiliary field. When dropping the contribution of 
gravitinos (tp^ = 0), there is the following formula for the 
superfield Lagrangian C(x, 9): 



S = J d 4 xd 2 9£C = J d i xe[BC 1 (x) + £ 2 {x)] , 



where 



Ci{x) = C\, C 2 (x) = V 2 C\ 



(3) 



(4) 



Here the vertical bars denote the leading field compo- 
nents of superfields. In particular one has 



B 



1 



pi 



3M pI 3Mj 
2 , 1 ^BB 



-(S + iP), 



(5) 



(6) 



where M p \ = 2.435 x 10 18 GeV is the reduced Planck 
mass. 

Applying the formula ([3]) to the action ([1]), we obtain 
the bosonic part of the supersymmctric Lagrangian of 
Eq. (P in the form [Io[ 

L = 3XF(X) + Qi? + AXX^j F'(X) + H.c. , (7) 

where X = B/(3AI p \). Here and in the next Sees. Mil and 
II VI we focus on the further reduced case in which the aux- 
iliary field B is real, i.e. X — X , by ignoring the complex 



nature of the supergravity superfields. In physical terms, 
it amounts to dropping the bosonic degree of freedom 
which is the pseudo-scalar supcrpartner of scalaron (see 
Sec. [Vj, in addition to dropping all fermionic degrees of 
freedom. Then the Lagrangian ([7]) reduces to 



L = 6XF{X) + 2 ( -R + AX 2 F'(X) . 



(8) 



Variation of the action §E§ with respect to X gives rise 
to the following algebraic relation between X and R: 

3F(X) + 11XF'(X) + Qi? + AX 2 ^j F"(X) = . (9) 

It is natural to expand the function F(1Z) into power 
series of 1Z, i.e. F(TZ) = J27=o Ci ^^ with some (generi- 
cally complex) coefficients Cj . Since we are not interested 
in the parity- violating terms in this paper, we assume 
that all the coefficients c, are real. For example, the 
choice F(TZ) = fo — fiTZ/2 gives rise to the standard 
Af = 1 supergravity with a negative cosmological con- 
stant. Slow- roll inflation can be realized by the following 
function 

F(TZ) = -\hTl + \f 2 n 2 - ±f 3 K 3 , (10) 

where /i 2,3 are positive constants having dimensions of 
[mass] 2 , [mass] 1 , and [mass] respectively. In this case 
the Lagrangian (JSj) yields 



L 



-\hR+\hBX- (jh + ^fsR 



+n.hx 3 - bf 3 x 4 . 

It follows from Eq. © that 



(11) 



-w 3 x2+ ^ R+R ^ x -w s R =°> ^ 



1 



h 



where 



Ro — 



21/i 

h 



(13) 



For the stability of the bosonic embedding given above 
we require the condition F'(X) < 0, which translates 
into f 3 X 2 - 2f 2 X + /i > [3. In order for this relation 
to hold for any real value of X, we have to require the 
condition 1 



/I < /1/3 



(14) 



For the successful inflation the second term on the 
r.h.s. of Eq. (|10j) needs to be suppressed relative to 



1 The stronger condition /| <S /1/3 was used in Ref. Il4l . 
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the third term. From the amplitude of the temperature 
anisotropies observed in CMB we require [l4| 

h > i • (is) 

In order to avoid large quantum corrections, the scalaron 
mass in the low-curvature regime [defined in Eq. 
needs to be smaller than the Planck mass, so that [L 



/I » h 



(16) 



There are two asymptotic regimes characterized by (A) 
\R\ > R and (B) \R\ < R . The former one is the high- 
curvature regime related to the generation of large-scale 
temperature anisotropies observed in CMB. The latter 
one is the low-curvature regime related to reheating after 
inflation. In the following we derive the approximate 
expression of A, as well as the Lagrangian L in those 
two different regimes, improving the results of Ref. [l4j [. 

. (A) \R\ » R 

In this case the first and third terms on the l.h.s. of 
Eq. (|12p correspond to the dominant contributions, 
so that the 0-th order solution satisfies Ag + (R + 
i?o)Ao/30 = 0. In order to connect this solution to 
the one in the regime \R\ <C Rq wc require Xq < 0, 
so that 



An = - 



R + Rq 



30 



(17) 



We regard the second and fourth terms in Eq. (TT2")) 
as the perturbations to the leading-order solution 
(TP?]). Setting A = A + 5X (\6X\ < |A |) and 
expanding Eq. ([12]) up to the linear order in 5X, 
the perturbation SX can be expressed in terms of 
/ 2 . This process leads to the following solution: 



A 



R + Rq 



30 



h 



13R + 33i? 



2 20f 3 (R + R ) - 33/ 2v /-30(i? + i? ) 



(18) 



Substituting this relation into Eq. (|11|) and expand- 
ing it at linear order in / 2 yields 



1 

180 



f3(R 2 + R 2 ) 



A 

900 



(9i? - URq)^^^ + R ) . (19) 



In the limit \R\ 3> Ro this Lagrangian reduces to 

30 



L ~ — -fiR+ — f 3 R 2 
RV 180 J 



100 



f-2(-R) 



3/2 



(20) 



. (B) \R\ « R Q 

In this regime the dominant contributions in 
Eq. (fT2"j) correspond to the third and fourth terms, 
so that 0-th order solution is given by 



An 



hR 



f 3 (R + R ) ' 



(21) 



which is negative. Taking into account the first and 
second terms as the perturbations to An, wc obtain 
the following solution: 



X 



hR 



MR + Ro) 



1 + ? 



f$R(13R + 33i? ) 



2 fl(R + R ) 3 - 9f$R(R + IIRq) 



.(22) 



Substituting this solution into Eq. (|11[) and expand- 
ing it in terms of / 2 up to the order of /| yields 



L 



If/? + ft R2 , ti(6R+llRa)R 3 
3 11 + 3/ 3 (i? + i?o) 



fi(R + R ) 4 



(23) 



Note that there is also the intermediate regime char- 
acterized by \R + Ro\ <^ Ro, i.e. when R is close to — Ro 
[Tl | . In this case the first and fourth terms in Eq. (fT2|) are 
the dominant contributions. Picking up the next-order 
contributions as well, we have 



X 



fM 1 ' 3 ^^. (24) 
V30/ 3 y 1 ' 20/ 3 V ' 



In this regime the Lagrangian (|11[) includes the term pro- 
portional to R alone. This transient era is followed by the 
reheating epoch characterized by the Lagrangian (|23l) . 

In order to recover the standard behaviour of Gen- 
eral Relativity in the low-energy regime we require that 
fi = 3A/pj/2. The mass squared of the scalar degree of 
freedom is given by m? = l/(3/"(i?)), where f(R) is rc- 

\R\ <§; Rq we have 



lated with L(R) as L(R) = -M* l f(R)/2. In the limit 



21/iMp 2 , _ 63M, 



pi 



4/| 8/| 



(25) 



In the high-curvature regime given by the Lagrangian 
(l20l) the scalaron mass squared is 



M 



15M| 

h 



(26) 



Then the constants /12 3 can be expressed by using the 
three mass scales M p i, to, and M, as 




15M p 2 j 



(27) 



The conditions dH|), (fT5|) . and ([!()]) translate into 



to > J — M , M <C M p i , to < M p i , (28) 
respectively. 
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III. INFLATIONARY DYNAMICS 

Let us study the dynamics of inflation for the f(R) 
model introduced in SecHH Here we are interested in the 
high-energy regime (A) satisfying the condition \R\ ^> 
Ro- 

We consider the flat Friedmann-Lemaitre-Robcrtson- 
Walkcr (FLRW) background described by the line ele- 
ment ds 2 = dt 2 — a 2 (t)dx 2 , where a{t) is the scale factor 
with the cosmic time t. The Ricci scalar is given by 
R = -6(2H 2 + H), where H = a/a is the Hubble pa- 
rameter (the dots stand for the derivatives with respect 
to t). 

In order to study the dynamics of inflation, it is conve- 
nient to introduce the following dimcnsionlcss functions: 



In the following we carry out the linear expansion in 
terms of the variables ti (i — 1,2,3), a, /3, and s = 
H/(HH). 

For the Lagrangian (|3U)) we have 



4H 2 ( 27V35 3 n 1 
8H 3 / 27^ 1 \ 



Then the variable C2 is given by 



<2 



/ , 27V35 3 . 1 1 

-ei 1 a a H — e-i -I — s 

1 I 800 40^ 2 1 4 



(35) 
(36) 

(37) 



M 2 
mH 



m 2 



(29) 



By using Eqs. (fL3")l and (|2"7|) . i? can be expressed as 
Ro = 21A/ 2 /10. During inflation the functions (|29|) 
should satisfy the conditions a <C 1 and ^ < 1 (sec 
below). In Eq. (Q2J} the term f 3 R 2 /180 is the domi- 
nant contribution during inflation. Hence, we neglect the 
higher-order terms beyond that of the first (linear) order 
in a and (3. Then the Lagrangian f(R) = -2L{R)/M^ 
following from (fll))) is given by 2 

3 R 2 3VW5(_R)^ 
^ 10 _ 6M 2 100 m~ " (30) 

We assume that the Lagrangian (f3U|) is valid by the end 
of inflation. 

In the flat FLRW spacetime the field equations of mo- 
tion are 3 



3TH 2 = (f -RF)/2-3HF, 
-2TH = f-HT, 



(31) 
(32) 



where T = f'(R). We introduce the slow- roll parameters 




H 
H 2 



* = e * = w> (33) 



which satisfy < 1 (i = 1,2,3). From Eq. it 
follows that 



£i = -£2(1 - £3) • 



(34) 



2 There are also some other approaches that give rise to the infla- 
tionary solution similar to that in the Starobinsky's f(R) model. 
Those include a large non-minimal coupling to the inflaton field 
[l6l | and a macroscopic theory of the quantum vacuum in terms 
of conserved relativistic charges [l7ll . 

3 Compared to the field equations given in Ref. 0] we only need 
to change R — v —R and / — v — /, while T is unchanged. 



Comparing this with Eq. (|34[) , we obtain 



£3 



27V35 3 „ 1 1 

a p H — ei H — s . 

800 40 M 2 1 4 



Similarly, Eq. (f3"Tj) gives the following relations: 

3V35 1 

ei = a H p, 

200 20 

and 



£2 



3V35 _ 1 
200 " 20^' 



Equation (|39[) is equivalent to 
Sv^M 2 



H 



200 m 



10m 



3V35 



(38) 



(39) 



(40) 



(41) 



This differential equation can be easily integrated. It 
yields 



H(t) = H t 



10771 



3V35 



exp 



f2 



3V35 M 
200 ~m 



iU-t) 



10m 



3V35 
(42) 

where Hi is the initial value of H at t = U . So we find 



(43) 
we obtain 



3V35 

s = a . 

200 

Substituting Eqs. §5§§ and (gSJ into Eq. 



£3 



3V35 1 
W a_ 20^' 



(44) 



The end of inflation (/j = ij) is identified by the condi- 
tion ei = 1. By using the solution (|42j) . we have 



t f, 20 ° TO lnf 63M 

1 ' ~~ 3%/35M 2 n V80m(3V35 J ff l + 10m) 

800 / m \ 2 



1 



SUU /my 
~63~ VAfJ 



1 



1600 /m\ 2 



63 VM 



.(45) 
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We define the number of e-foldings from the onset of 
inflation (t = ij) to the end of inflation (t = tf) as 

N(U) = f* f Hdt. From Eqs. (g2]) and (gSJ) we can ex- 
press N(ti) in terms of Hi, M, and m. The number of 
c-foldings N corresponding to the time t can be derived 
by replacing Hi in the expression of N(ti) for H. It fol- 
lows that 



N = 



1 



3a(80V35 - 21a - V7(63a 2 + 1600/3)) 



126a 2 _ 
400/3(81n2 + 31n5) + 800/3 



In 



\/7(63a 2 + 800/3) + 21aV63a 2 + 1600/3 



21a + 2V35/3 



(46) 



In the limit a -> one has N -> 10//3 - 1/2, i.e. 
/3 -> 20/(2iV + 1). In this case the i? 2 /(6M 2 ) term 
in the Lagrangian (|30[) dominates over the dynamics of 
inflation, which corresponds to the Starobinsky's /(i?) 
model [l|. In another limit j3 — > it follows that 
N -> 40^35/(210) - 1, i.e. a -> 40V35/[21(A r + 1)]. 
Then we obtain the following bounds on a and /3: 



< a < 



40V35 
21(iV + 1) 



< p < 



20 



2iV + 1 



(47) 



In order to realize inflation with TV = 60, for example, 
the two variables need to be in the range < a < 0.185 
and < j3 < 0.165. For the number of e-foldings relevant 
to the CMB temperature anisotropies (50 < N < 60) the 
slow- roll parameters given in Eqs. ([3"9"]l . (|4"0"|). (|H} are 
much smaller than unity, so that the slow-roll approxi- 
mation employed above is justified. 



IV. OBSERVATIONAL TESTS 

Let us study whether the f(R) model ([3"U]) can sat- 
isfy observational constraints of the CMB temperature 
anisotropies. The power spectra of scalar and tensor per- 
turbations generated during inflation based on f(R) the- 
ories have been calculated in Ref. 01 (see also Ref. 0). 

The scalar power spectrum of the curvature perturba- 
tion is given by 



— ( — 



Using Eqs. ([35]), ([39]), and ([i0]). it follows that 



1250 ( M\ ( i— \- 2 



(48) 



(49) 



where, in the expression of J- , we have neglected the 
terms a and (3 relative to 1. Using the WMAP nor- 
malization V s = 2.4 x 10 -9 at the pivot wave number 
fco = 0.002 Mpc -1 [lH, the mass M is constrained to be 



M ~ 7.5 x 10" 6 (3V35a + 10/3 j M pl . 



(50) 



0.0010 




0.01 

a 



Figure 1: The two masses M and m versus the variable a in 
the regime 10 -4 < a < 0.18 for the number of e-foldings N = 
55. We also show the upper bound a max = 0.201 determined 
by Eq. (|47p . M is weakly dependent on a with the order of 
10 Mpi, whereas m strongly depends on a. The condition 
m > a/7/20 M is satisfied for a < 0.178. 



In the limit that a -> and /3 20/(2iV + 1) we 
have M/Mpi = 7.5 x 10~ 4 /(7V + 1/2). In another limit 
a -> 40V^5/[21(A r + 1)] and /3 -> it follows that 
M/M p i = 1.5 x 10~ 3 /(/V + l). In the intermediate regime 
characterized by (|47p we can numerically find the values 
of a and /3 for given 2V satisfying the constraint (|46p . 
which allows us to evaluate M from Eq. ([50]) . From 
Eq. ([2^]) the mass scale m is also known by the relation 
m = (y/]5/a)M. 

In Fig. [1] wc plot M and m versus a in the regime 
10~ 4 < a < 0.18 for /V = 55. In this case a is bounded 
to be < a < 0.201 from Eq. (|i7]). The mass M weakly 
depends on a with the order of 10 _5 M p i, while m changes 
significantly depending on the values of a. For a much 
smaller than 1 we have m 3> M, while m is of the same 
order as M for a > 0.1. Wc recall that there is the 
condition m > ^7/20 M coming from the requirement 
(ITU . For /V = 55 this condition gives the upper bound 
a < 0.178. 

The scalar spectral index n s is defined by n s = 1 + 
dlnVs/dhik, which is evaluated at the Hubble radius 
crossing k = aH (where k is a comoving wave number) 
HIUl. In f(R) gravity this is given by 



1 - 4ei + 2e 2 - 2e 3 . 



On using Eqs. (p39l. ([40]) and (EH]) , we obtain 



:/3- 



3V35 
100 



(51) 



(52) 
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The tensor power spectrum is given by Q 



V, = 



H 



TT J~ \M pl 

From Eqs. (|48f and (|53|) the tensor-to-scalar ratio is 
_ Vt ,„ 2 3 



(53) 



2500 



3V35a + 10/3 



(54) 



In the limit a -> and /3 -> 20/(27V + 1) the observ- 
ables dsn and (El) reduce to 



n s (a -> 0) 
r(a 0) 



1 



27V- 
48 



1 ' 



(27V + l) 2 ' 



(55) 
(56) 



which agree with those in the Starobinsky's f(R) model 
@. For TV = 55 one has n s (a — > 0) — 0.964 and r(a -> 

0) = 3.896 x 10~ 3 . In another limit a -> 40 v / 35/[21(7V + 

1) ] and j3 — > it follows that 



n s (/3^0) = 1- 
r(P -> 0) = 



7V + 1 ' 



(TV + 1) 2 



(57) 
(58) 



For TV = 55 one has n s (/3 -> 0) = 0.964 and r(/3 -> 
0) = 1.531 x 10~ 2 . While the scalar spectral indices 
(|55|) and (|57|) are practically identical for TV ^ 1, r(/3 — >■ 

0) is about four times as large as r(a — > 0). For the 
intermediate values of a between and 40-\/35/[21(7V + 

1) ] we need to numerically derive /3 satisfying Eq. (|46[) 
for given TV. Then n s and r are known from Eqs. (|52[) 
and |(51J). 

In Fig. [2] we plot the theoretical values of n s and r in 
the (n s , r) plane for TV = 50, 60, 70 together with the la 
and 2(7 observational contours constrained by the joint 
data analysis of WMAP7 fl8| . Baryon Acoustic Oscilla- 
tions (BAO) [2l|, and the Hubble constant measurement 
(HST) [22j] ■ The observational bounds are derive d by 
using the standard consistency relation r = — 8n t [20|, 
where n t = d In Vt / d In k is the tensor spectral index. In 
f(R) gravity this relation also holds by using the equiva- 
lence of the power spectra between the Jordan and Ein- 
stein frames Q (see also Refs. [23j). 

The Starobinsky's f(R) model, which corresponds to 
the limit a — > with the observables given in Eqs. (J55J 
and (|56p . is well within the current observational bound. 
In the regime a < (3 one has /3 ~ 20/(27V+ 1) - V35a/5, 
so that 



n s (a < /?) 
r(a < 



1 



35 



2TV+1 100 



a , 



(27V +1) 2 



35 (27V+ 1) 
^00 ' 



(59) 



.(60) 



This shows that both n s and r increase for larger a sat- 
isfying the condition a <C /?■ As we see in Fig. [2J ?i s 




0.92 0.94 0.96 0.98 

n 



Figure 2: The three thick lines show the theoretical values of 
n s and r for TV = 50, 60, 70 with a ranging in the region (|47[) . 
The thin solid curves are the la (inside) and 2a (outside) 
observational contours constrained by the joint data analysis 
of WMAP7, BAO, and HST. For a -> 0, n s and r are given 
by Eqs. (|55p and (|56p . In the limit j3 — >• 0, n a and r approach 
the values given in Eqs. (|57[) and (I58|l . 



switches to decrease at some value of a, whereas r con- 
tinuously grows toward the asymptotic value given in 
Eq. ((53). 

From Fig. d we find that the f(R) model (0 in 
which a is in the range (|47[) is inside the Ict observa- 
tional contour. The condition m > v/7/20M provides 
the constraints a < 0.194, a < 0.165, a < 0.143 for 
TV = 50, 60, 70 respectively, while the bound (|47|) in each 
case corresponds to a < 0.221, a < 0.185, a < 0.159. 
When TV = 60 the scalar spectral index and the tensor- 
to-scalar ratio are n s = 0.969, r = 0.0110 for a = 0.165 
and n s = 0.967, r = 0.0129 for a = 0.185, which are not 
very different from each other. For the background in 
which inflation is sustained with the number of e-foldings 
N > 50 the model is consistent with the current obser- 
vations. 

Note that the nonlinear parameter /nl of the scalar 
non-Gaussianities is of the order of the slow-roll parame- 
ters in f(R) gravity [24| . Hence, in current observations, 
this docs not provide additional constraints to those stud- 
ied above. 



V. EFFECTIVE SCALAR POTENTIAL 

In this section we derive the effective scalar potential 
and the kinetic terms of a complex scalaron field in the 
low-energy regime (B) characterized by \R\ <C Ro- In 
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doing so, let us return to the original F(TZ) supergravity 
action (TP and perform the superiield Legendre transfor- 
mation |25| ■ We temporarily set M p i = 1 to simplify our 
calculations. It yields the equivalent action 



S 



d 4 xd 2 6E 



-yn + z{y)} + u.c, (6i) 



where we have introduced the new covariantly chiral su- 
pcrficld y and the new holomorphic function Z(y) re- 
lated to the function F as 



F(ii) = -ny{n) + z(y(n)) . 



(62) 



The equation of motion of the supcrficld y, which follows 
from the variation of the action (|6Tj) with respect to y, 
has the algebraic form 



K = z\y) 



(63) 



so that the function y(lZ) is obtained by inverting the 
function Z' . Substituting the solution y(7V) back into 
the action (|61[) yields the original action ([1]) because of 
Eq. 162]). We also find 



y = -F'(n) . 



(64) 



The inverse function 7Z(y) always exists under the physi- 
cal condition F'(JZ) ^ 0. As regards the F- function (|T0l) . 
Eq. (jM)) yields a quadratic equation with respect to 1Z, 
whose solution is 



n{y) 



UM 2 
20m 



80m 2 
21Af 2 



(y - 3/4) 



(65) 



where we have used the paramctrization (f2"T|). Equation 
(p5|) is also valid for the leading complex scalar field com- 
ponents 1Z\ = B/3 = X and y\ = Y, where Y is the 
complex scalaron field. 

The kinetic terms of y are obtained by using the iden- 
tity 



d i xd 2 8£ yK + YL.c. = d A xd A 9E- 1 {y + y), (66) 



where E 1 is the full curved superspace density [l| 
Therefore, the Kahler potential reads 



k = -3 in (y + y) 

It gives rise to the kinetic terms 
d 2 K 



(67) 



£kin — 



= 3 



oyay 



d u Yd^Y 



d^Yd^Y 



y=v 



[Y + Yf 



4y2 



(68) 



where we have used the notation Y = y + iz in terms 
of the two real fields y and z. The imaginary com- 
ponent z corresponds to a pseudo-scalar field. The ki- 
netic terms (|68|) represent the non-linear sigma model 



|26| with the hyperbolic target space of (real) dimension 
two, whose metric is known as the standard Poincare 
metric. The kinetic terms are invariant under arbitrary 
rescalings Y — > AY with constant parameter A ^ 0. 

The effective scalar potential V(Y, Y) of a complex 
scalaron Y in the regime (B), where supergravity decou- 
ples (it corresponds to rigid supersymmetry) , is easily de- 
rived from Eq. (|61|) when keeping only scalars (i.e. ignor- 
ing their spacetime derivatives together with all fcrmionic 
contributions) and eliminating the auxiliary fields, near 
the minimum of the scalar potential. We find 



v/ = ||z'(r)| 2 = ||^(y)| 2 , 

which gives rise to the chiral superpotential 



w(y) = \i Y z (y) 



(69) 



(70) 



The superfield equations (|67|) and (|70|) are model- 
independent, i.e. they apply to any function F(7Vj in 
the large M p \ limit, near the minimum of the scalar po- 
tential with the vanishing cosmological constant. 

There is no field redefinition that would bring all the 
kinetic terms (|68p to the free form. The canonical (free) 
kinetic term of a real scalaron y alone can be obtained 
via the field redefinition 



y = J 4cxp(- v /273 0) ■ 



(71) 



The scalaron potential vanishes at y = 3/4. Demand- 
ing that this minimum corresponds to <j> = 0, we have 
A = 3/4 and hence y = (3/4) exp(— \/2/3 <f>). Defining a 
rescaled field \ as X — \/8/3 z, the kinetic term (|68)l can 
be written as 



kin 



(72) 



Here and in what follows we restore the reduced Planck 
mass M p \. 

The total potential (|69|) including both the fields (j> and 
X is given by 



V(cf>,x) = 
where 

B{4>) 

c(x) 



U7M 4 Ml l 
400m 2 



1 



20m 2 



7M 2 
)m 2 [3 x 



^B{4>) + iC( X )-l 



-y/2/3 0/M p i 



21Af 2 



(73) 

(74) 
(75) 



In order to express (|T3")) in a more convenient form we 
write ^JB{4>) + iC(x) — P + iq, where p and q are real. 
This gives the relations p 2 — q 2 = B{<f>) and 2pq = C(x)- 
Solving these equations for p, we find 



1 



B{<t>) + VWW+c^x) 



1/2 



(76) 
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w here w c have chosen the solution p > to recover p = 
y/B{<f>) for B{<f>) > in the limit C(x) -> 0. Then the 
field potential (f73|) reads 



V{cj>, X ) 



147M 4 M 



pi 



400m2 WW WW 



1/2 



(77) 



In the absence of the pseudo-scalar x the potential ([77] 
reduces to 



V(</>) 



147M 4 M 2 1 
400m 2 



1 + |S(0)| -V2{S(0) + |S(0)|} 1/2 

(78) 

For the field <fi satisfying the condition B(<fi) < it follows 
that 



V{4>) 



21 
20 



M 2 A/ 2 (l - e -y/m*/M#j , (79) 



which approaches the constant V(<j>) — > 21M 2 M 2 1 /20 in 
the limit (f> — > oo. Defining the slow- roll parameter ey = 
(M2/2)(l/ /F) 2 , we have 



3(1 -x) 2 ' 



/2/3 <j>/M p 



(80) 



The end of inflation is characterized by the criterion ey = 
1. This gives x f = e-Vystf/M,* = (3-^/2 and hence 
(f) f = 0.558M p i. For m > M the condition #(</>/) < is 
satisfied, so that the potential ([79"]) is valid at the end of 
inflation. If m is close to the border value ^/7/20M - 
see Eq. (|28p — , then the potential ((79")) is already invalid 
at the end of inflation. 

For small <fr satisfying the condition B{4>) > the po- 
tential (|78p reads 



V(4>) 



147M 4 M 2 1 
400m 2 



1 



/ 20m 2 / 



-y/2/34/Mpi _ 1 



(81) 

In this case Taylor expansion around <j> = gives rise to 
the leading-order contribution V(<f>) = m 2 2 /2. Reheat- 
ing occurs around the potential minimum through the 
oscillation of the canonical field <j>. 

The total effective potential involving the interaction 
between the fields <j> and x is given by Eq. ([77]) . Expand- 
ing the potential (|77|) around <j) = x = and picking up 
the terms up to fourth-order in the fields, wc obtain 



V(<j>,x) 



2 j.2 



-m 4> 



% /6m 2 (10m 2 - 7M 2 ) 



(1500m 4 



5\/6m 4 
"21M 2 M p i 



<f>X 



42M 2 M p i 
2 5m 4 (10m 2 



147M 4 M 2 ! 



™ 2 W 



1260m 2 M 2 + 343M 4 )m 2 4 1 2 2 25m 6 4 



1764M 4 M 2 ! 



+ 2 m X 49M 4 M 2 



pi 



(82) 



The scalaron </> is coupled to to the pseudo-scalar x 
through the interaction given in the second line of 
Eq. 



VI. PREHEATING AFTER INFLATION 

We study the dynamics of preheating for the two-field 
system described by the kinetic term (|72|) and the effec- 
tive potential ([77]) . The background equations of motion 
on the flat FLRW background are 

3M^H 2 = <j> 2 /2 + e 2b X 2 /2 + V, (83) 
+ 3H<j> + V<p- b^e 2b X 2 = , (84) 
X+(3H + 2b^) X + e- 2b V x =0, (85) 

where b{4>) = \/2/3 4>/M p i and " ^" represents a partial 
derivative with respect to </>. 

In Fourier space the field perturbations 8<f) k and $Xk 
with the comoving wave number k obey the following 



equations (see Refs. [27J for related works): 

5<j> k + 3HS<p k + [k 2 /a 2 + V M - {2b\ + b, H )e 2b X 2 }Hk 
= -V^Sxk + 2b^e 2b xS Xk , (86) 
5 X k + (3H + 2b^)S Xk + (k 2 /a 2 + e~ 2b V, xx )5xk 
= -e- 2b {V^ x - 2b.^V x + 2b, u e 2b 4>x)Hk - 2&,**ty fc . 

(87) 

The derivative V x of the potential ([77]) vanishes at 
X = ±Xc, where 



Xc = 



/210M 
20m 



1 _ e -\/2/3 <t>/M pl 



21 / M 



80 V m 



n 1/2 



M, 



The local minima exist in the x direction provided that 

-n -1 



> \ -In 
V 2 



1 



21 M 



80 V m 



M p i = C , (89) 
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Figure 3: Effective potential (JTZj) for m = 1.14 x 10" 4 M p i 
and M = 1.62 x 10~ 5 M p i in the region -0.1 < <f>/M p i < 0.1 
and -0.1 < x/Mpi < 0.1. 



X — Xc- After </> drops below fa., the field x approaches 
the global minimum at % = 0. In Fig. U we show one 
example for the evolution of the background fields <j> and 
X with the same values of m and M as those in Fig. [3] 
The energy density of the field x catches up to that of 
the infiaton around the onset of reheating. 

As we see in Eq. (|89")l . the critical field value <fi c gets 
smaller for increasing m/M. Hence, for larger m/M , the 
potential (|77|) possesses the local minima at x = ±Xc for 
a wider range of 0. The potential in the region \x\ < Xc 
can be flat enough to lead to inflation by the slow-roll 
evolution of the field Xj even if <p is smaller than fat = 
0.558Mpi. For larger ratio m/M inflation ends with the 
field value much smaller than fat. If m/M = 20 and 
m/M = 83, for example, the amplitudes of the field 4> 
at the onset of oscillations are fa = 1.5 x 10 _2 M p i and 
fa = 5.0 x 10 _3 M p i, respectively. 

Let us consider the regime where the condition 



/ my 

\m) 



101 



M, 



< 1 



i>i 



(90) 




30 



mt 



Figure 4: Evolution of the background fields <jr and x (both 
are normalized by Mp\) for m = 1.14 x 10 _4 Af p i and M = 
1.62 x 10 _5 M p i with the initial conditions </> = 0.55M p i, x = 
10" 3 Af p i, 4> = -1-6 x 10" 2 mM p i, and x = 1.5 x 10" 3 mM p i. 



whereas they disappear for <p < fa. In Fig. |3]we plot the 
potential (|77|) with respect to <j> and x f° r 171 = 1-14 x 
10" 4 Af pl and M = 1.62 x 10" 5 M pl . Since fa = 6.5 x 
10 _3 M p i in this case, the potential has the local minima 
in the x direction for <j> > 6.5 x 10~ 3 M p i. From Eq. ([88]) 
the field value Xc increases for larger fa For the model 
parameters used in Fig. [3J for example, one has Xc = 
0.028M p i at <j> = 0.1M p i and Xc = 0.059M pl at = 
0.5M p i. 

If the initial conditions of the fields are < X < Xc and 
4> > 4> c , the field x grows toward the local minimum at 



is satisfied. Then the potential (|82|) is approximately 
given by V(fa x) — m 2 (j> 2 /2 + m 2 x 2 /2, in which case 
both <p and x have the same mass m. This gives rise to 
the matter-dominated epoch (where H = 2/(3£)) driven 
by the oscillations of two massive scalar fields. From 
Eq. dBU we have that $ + (2/i)0 + m 2 4> ~ 0, whose 
solution is 



0W - ^ — T^i sin(TTit) . 
zmt 



(91) 



Here the initial field value fa corresponds to the time 
U = 7r /(2m). 

In order to discuss the dynamics of the field perturba- 
tions in Eqs. ([5B"f and ([57)) we define the two frequencies 
cj0 and uj x , as = k 2 /a 2 + - (2b 2 (jl + b 4 



b)e 2h x 2 and 



= k 2 /a 2 + e 



-2b 



V 



As long as the condition (f50"| is 



satisfied, it is sufficient to pick up the terms up to cubic 
order in fields. It then follows that 



w,7 



a 2 
k 2 



v / 6m 2 (10m 2 - 7M 2 ) 



pl 



+ m 2 e- 2b 



7M 2 M ] 

lOV6m i 
21A'/ 2 M pl ' 



-26 



4>. 



(92) 
(93) 



where, in Eq 
the term — (26 2 



we have neglected the contribution of 

b^)e 2 



" 2b X 2 



We introduce the rcscalcd fields 5(fk 



,3/2; 



and 



SXk = a 3 / 2 e b Sxk to estimate the growth of perturba- 
tions in the regime (|90"]) . Neglecting the contributions of 
the r.h.s. of Eqs. ([86]) and (|87| and also using the ap- 
proximation e _2b ~ 1 in the regime H <C m, the field 
perturbations Scpk and SXk obey the following equations 



j- 2 6 V k + [A k 



, cos(2z)] 5 Vk ^ , (94) 



■^SXk + [A k - 2q x cos(2z)] SX k s , (95) 
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Figure 5: Evolution of the field perturbations 8(f> k = 
k 3 / 2 5(f)k/M p i and S\k = k 3 ^ 2 8xk/M p i with the wave num- 
ber k — m for m = 1.16 x 10~ 3 Af p i and M = 1.39 x 10" 5 M p i. 
We choose the background initial conditions <j> = 0.1M p i, 
X = 1.0 x 10" 3 Af p i, 4> = -8.48 x 10" 4 mM p i, and x 

i . : x x icr 5 mM p i. 




10' 6 I I 
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mt 

Figure 6: Evolution of the field perturbations with the wave 
number k = m for m = 2.89 x 10" 4 Af p i and M = 1.46 x 
10~ M p i. We choose the background initial conditions (j> = 
O.IA'V, x = 1.0 x 10" 3 Af p i, <t> = -7.35 x 10" 3 mA/ p i, and 
X = 6.85 x 10~ 4 mA/ pl . 



where 2z = mt + w/2. The quantities Ak, q^,, and q x are 
given by 



20V6 / 7M 2 \ /mx2 & n/2 

ff * = "T" 1 1_ io^JImJ Iw' (97) 

20V6^ 2j W2 ; (Q8) 
yx 21 VM/ M p i mt ' v y 

which are time-dependent. 

Equations (|94| and (|95|) are the so-called Mathieu 
equations describing parametric resonance caused b y th e 
oscillation of the field (f> [28|,[2!| (see also the review pfj|). 
In the regime (|90|) both and q x are smaller than 1 for 
t > ti = 7T /(2m). In this case the resonance occurs in 
narrow bands near Au = Z 2 , where I = 1, 2, • • • [29l. [30j . 
As the physical momentum fc/a redshifts away, the field 
perturbations approach the instability band at Ak = 4. 
Although 6(fik and SXk can be amplified for Ak — 4 and 
50 ~ li <?x ~ 1' thi s narrow parametric resonance is not 
efficient enough to lead to the growth of 5(f>k and (>Xfc 
against the Hubble friction [29[. 

If the initial field <fii satisfies the condition 
(m/M) 2 \4>i\/M v \ ^> 1, the quantities and q x are 
much larger than 1 at the onset of reheating. This 
corresponds to the so-called broad resonance regime 
[29j in which the perturbations 5<pk and S\k can grow 
even against the Hubble friction. We caution, however, 
that Eqs. (|M|) and arc no longer valid because the 



background solution (jUTj) is subject to change due to the 
effect of higher-order terms in the potential ([77]) . Still, 
the non-adiabatic particle production occurs around 
the potential minimum ((/> = 0) [2ll . In this region the 
dominant contribution to the potential is the quadratic 
term m 2 (f) 2 /2. Hence it is expected that preheating can 
be efficient for the values of q^ and q x much larger than 
1 at the onset of the field oscillations. 

Numerically we solve the perturbations equations ([55)1 
and (|87p together with the background equations (|83p . 
([55)1. and (@SJ for the full potential (J77J) without using 
the approximate expression (|82|) . In Figs. [5] and [6] we 
plot the evolution of the field perturbations 64>k and Sxk 
with the wave number k = m for two different choices of 
the parameters m and M (which are constrained by the 
WM AP normalization in Fig. [1} . The initial conditions of 
the perturbations are chosen to recover the vacuum state 
characterized by SifkiU) = e _lw * ti /^/2o;0 and SXk(ti) — 

Figure [S] corresponds to the mass scales m = 1.16 x 
10" 3 M pl and M = 1.39 x 10" 5 M pl , i.e., the ratio m/M = 
83. The field value at the onset of oscillations is found 
to be 0j = 5.0 x 10 -3 Af p i, in which case q^>{U) = 244 
and q x (ti) = 81. Figure [S] shows that both 54>k and 8xk 
grow rapidly by the broad parametric resonance. The 
growth of the field perturbations ends when and q x 
drop below 1. 

Figure [6] corresponds to the ratio m/M = 20, in which 
case <j)i = 1.5 x 10 _2 M pl , q$(U) = 41, and q x (ti) = 4.6. 
Compared to the evolution in Fig. (SJ preheating is less 
efficient because of the smaller values of q^iti) and q x (ti). 
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The parameter to control the efficiency of preheating is 
the mass ratio m/M. For larger m/M the creation of 
particles tends to be more significant. For the mass m 
smaller than 10~ M p \ the field perturbations 6<fik and 5\k 
hardly grow against the Hubble friction because they are 
not in the broad resonance regime. 

In our numerical simulations we did not take into ac- 
count the rescattering effect between different modes of 
the particles. The lattice simulation [3J HH is required 
to deal with this problem. It will be of interest to see 
how the non-linear effect can affect the evolution of per- 
turbations at the final stage of preheating. 

VII. CONCLUSIONS 

We have studied the viability of the f(R) inflationary 
scenario in the context of F(1Z) supergravity. In the high- 
energy regime characterized by the condition \R\ Rq 
there is a correction of the form (— R) 3 / 2 /m to the func- 
tion f(R) = 3R/10- R 2 /(6M 2 ). Introducing the dimen- 
sionless functions a and (3 in Eqs. (|29|). we showed that 
these are constrained to be in the range (|T7)) to realize 
inflation with the number of e-foldings N. 

The masses of the scalaron field in the regimes \R\ 
Rq and \R\ <C Rq are approximately given by M and 
m, respectively. From the WMAP normalization of the 
CMB temperature anisotropics we derived M and m as 
a function of a in Fig. Q] The weak dependence of M 
with respect to a means that the term —R 2 / (6M 2 ) needs 
to dominate over the correction (— R) 3 / 2 /m during in- 
flation. We also showed that the model is within the 
la observational contour constrained from the joint data 
analysis of WMAP 7, BAO, and HST, by evaluating the 
scalar spectral index n s and the tensor-to-scalar ratio r. 

In the presence of the pseudo-scalar field x coupled to 
the scalaron field (f> we derived the effective potential ([77)1 
and their kinetic energies (|72p in the low-energy regime 
(\R\ Rq). Provided that the condition (|59"|) is satisfied, 
the effective potential has two local minima at x = iXc- 
Around the global minimum at <fi = \ = the system is 
described by two massive scalar fields with other interac- 
tion terms given in Eq. ([52)1 . Even if \ is initially close 
to 0, x typically catches up to 4> around the onset of the 
field oscillations (see Fig.H]). 

In the regime where the field 4> is in the range (|M|) we 
showed that both the field perturbations 8<pk = a 3 ^ 2 S(f>k 
and SXk = a 3 / 2 e b 8xk obey the Mathieu equations ([94]) 
and (|9"5j) . This corresponds to the narrow resonance 
regime in which and q x are smaller than the order 
of unity. The broad resonance regime is characterized by 
the condition (m/M) 2 \<p\/M p \ 3> 1, but in this case the 
expansion (|82p of the effective potential around the min- 



imum is no longer valid. In order to confirm the presence 
of the broad resonance we numerically solved the pertur- 
bation equations ([55} and (|57|) for the full potential ([77]) . 
Indeed we found that preheating of both the perturba- 
tions 6cj>k and 8\k is efficient in this regime. As we see 
in Figs. [5] and [BJ the broad parametric resonance is more 
significant for larger values of m/M . 

Our results lend compelling support to the phenomeno- 
logical viability of the bosonic sector of F(1Z) super- 
gravity, in addition to its formal consistency. It is also 
worthwhile to recall that supergravity unifies bosons and 
fcrmions with General Relativity, highly constrains parti- 
cle spectrum and interactions, has the ideal candidate for 
a dark matter particle such as the lightest super-particle 
[33j], and can be deduced from quantum gravity such as 
superstring theory. The F(TZ) supergravity action ([1]) is 
truly chiral in superspace, so that it is expected to be 
protected against quantum corrections [l2| , which is im- 
portant for stabilizing the masses M and m in quantum 
theory. 

After the broad parametric resonance and the subse- 
quent rescattering at the final stage of preheating, decay 
of super-scalaron produces new particles including the 
visible sector. It can be studied perturbatively, when 
adding a supcrsymmetric matter action to the F(1Z) su- 
pergravity action in E q. CD ), along the standard lines, see 
e.g., Section 13 of Ref. [12| for a review. As is well known, 
the Starobinsky f(R) gravity model has the universal re- 
heating mechanism due to the coupling of scalaron to 
all matter [l], Q . The same applies to super-scalaron in 
matter coupled F(1Z) supergravity fl4j . 

Of course, our analysis of F{TZ) supergravity is still in- 
complete since it does not include fermions such as grav- 
itinos and inflatinos (gravitino is the fcrmionic supcrpart- 
ner of graviton, and inflatino is the fcrmionic supcrpart- 
ner of scalaron). We did not study in this paper further 
particle production to complete the reheating process. 
A calculation of decay rates and particle abundances af- 
ter preheating requires an extension of the action ([1} by 
some hidden sector to be responsible for supersymmetry 
breaking. It should be the subject of a separate investi- 
gation. 
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